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Inspired by recent experiments on 3 He films between one and two atoms thick, we consider a bilayer 
Hubbard model on a triangular lattice. Our results are obtained in the framework of a cluster 
dynamical mean- field calculation with a quantum Monte Carlo impurity solver. For appropriate 
model parameters, we observe an enhancement of the effective mass as the first layer approaches 
integer filling and the second remains only partially filled. At finite temperatures, this increase of 
the effective mass — or, equivalent ly, the decrease of the coherence temperature — leads to a crossover 
to a state where the first layer fermions localize, drop out of the Luttinger volume, and generate 
essentially free local moments. This finite temperature behavior is shown to be robust against 
the cluster size above some critical temperature. The zero temperature phase diagram, however, 
depends on the cluster topology. In particular, for clusters with an even number of unit cells, the 
growth of the effective mass is cut off by a first-order, orbital-selective Mott transition. 



I. INTRODUCTION 

The solidification of 3 He monolayers^ has been inter- 
preted as a density-driven Mott transition in which the 
effective mass diverges. 2 3 Below the critical density, the 
system is a metallic, nearly localized Fermi liquid; beyond 
the critical density, it is a solid, the magnetic properties 
of which are dominated by antiferromagnetic two-body 
exchange pr o cesses P"^ It is now possible to realize bilay- 
ers of 3 He (atop a frozen 4 He substrate, itself adsorbed 
onto graphite) with the special property that the sec- 
ond layer begins to form before the first has solidified. 7 
Since the first layer is close to a Mott transition, the 
3 He fermions in this layer are slow (i.e. heavy), whereas 
those in the second layer are fast. This combination of 
fast and slow dynamics — corresponding to wide and nar- 
row fermion conduction bands — is completely analogous 
to the situation in electronic heavy fermion materials, al- 
beit without the complication of crystal field and spin 
orbit effects. 

According to this picture, one expects, prior to solidifi- 
cation of the first layer, an enhanced effective mass and a 
Luttinger volume that counts both the first- and second- 
layer populations. Moreover, one naively anticipates that 
further 3 He deposition will eventually cause the effective 
mass to diverge, in coincidence with the solidification of 
the first layer. This solidification of the first layer can 
be interpreted either as an orbital-selective Mott transi- 
tion or, in the terminology of Kondo physics, as a Kondo 
breakdown in which the heavy particles drop out of the 
Luttinger volume. In experiment, the effective mass is 
indeed observed to increase as a function of the total 
3 He concentration, but its growth is interrupted by an 
intervening phase P The fact that this phase is ferromag- 
netic indicates that three-body exchange processes come 
to dominate in the solid phase of the first layer 

The motivation of this article is to consider a simple 
lattice model that goes a good way towards reproducing 
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FIG. 1. (a),(b) Stacking of billiard balls modeling of bi- 
layer 3 He, top and side view, with the 4 He substrate shown 
in white, (c) Tight-binding modeling with hoppings t c , £/, 
and V. (d) The hexagonal Brillouin zone of the triangular 
lattice, (e) The set of supercells, each consisting of N u unit 
cells, considered in this work. 



the essential features of the above experimental situation. 
As shown in Figs. [TJa) and [TJb), we adopt a stacking 
of billiard balls modeling of bilayer 3 He on a triangular 
lattice defined by a x = (1/2,^3/2,0) and a 2 = (1,0,0). 
Each unit cell accounts for two 3 He positions, rj = 
and r c = |ai — |a 2 + (0, 0, as), measured relative to the 
lattice. This geometry presupposes a particular stacking 
arrangement for the second 3 He layer. 

Our model can be viewed as a honeycomb lattice whose 
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inequivalent sites (corresponding to 3 He positions in the 
upper and lower layers) are populated by two species of 
fermion, which we label c and /. The tight-binding pa- 
rameters include a nearest-neighbor (interlayer) hopping 
V and next-nearest-neighbor (intralayer) hoppings t c and 
tf. See Fig. [IJc). With the inclusion of onsite Coulomb 
repulsion terms, the Hamiltonian reads 

f M/^c(k)-/i V(k) \fc K a\ 
U -2^ {%* J*,*) { V (k) e f (k)-fi) {f Ka J 

k,cr x / \ / 

+ C/ c ^(n c , i -l) 2 + [/ / ^(n / , i -l) 2 

(1) 

Here, the mixing element V(k) = V(3 + 27k) 1 / 2 and the 
dispersion e c (k) = —2t c ~f\r_ + are expressed in terms 
of the connection 7k = cos(k • ai) + cos(k • a 2 ) + cos[k • 
(a2 — ai)] of the underlying Bravais lattice. The operator 
^c,i = Xlcr c \ a c hv * s ^ ne local 3 He density in the upper 
layer. Similar definitions hold for £/(k) and hf^. 

Except for the complication of the layer stacking (and 
the resulting k-dependent hybridization), this bilayer 
Hubbard model reduces to the Periodic Anderson Model 
as t/ — >• 0, a limit in which the bare mass of the / 
fermions diverges. Similar models have been considered 
for the description of bilayer 3 He in Refs. [8] and |9] within 
a slave boson mean-field calculation. Here we go a signifi- 
cant step further and perform calculations within the cel- 
lular dynamical mean field theory (CDMFT)P^ approxi- 
mation (Sec.|H|. Our strategy is to systematically inves- 
tigate the model of Eq. Q as a function of the cluster 
size. In Sec. |III[ we will see that clusters with an odd 
number of unit cells have a radically different low-energy 
behavior than those with an even number. Given this 
situation, the extrapolation to the large cluster size limit 
is delicate and is relegated to the conclusions in Sec. |IVJ 
Part of this work has already appeared in a preprint 



II. CELLULAR DYNAMICAL MEAN FIELD 
THEORY 

By construction, the CDMFT approach exactly ac- 
counts for the temporal fluctuations at each site and 
thereby captures the physics of the local moments — both 
their formation and their screening via the Kondo effect. 
But the spatial fluctuations extend only over the simu- 
lation cluster; insofar as the true correlation length scale 
exceeds the linear size of the cluster, the results will suffer 
from finite size effects. To mitigate this, we have consid- 
ered various cluster topologies ranging from one unit cell 
(a single c and / site) to four unit cells as defined in 
Fig.^e). For a given supercell, the resulting single par- 
ticle Green function, G(K,zcj m ), is a 2N U x 2N U matrix 
with crystal momentum K in the Brillouin zone of the 
supercell lattice. The CDMFT calculation involves ne- 
glecting momentum conservation and thereby obtaining 
a K-independent self-energy S(zcj m ). This quantity is 



extracted from a cluster of N u unit cells embedded in a 
dynamical mean field that is determined self-consistently. 
We have solved this cluster problem using a standard 
Hirsch-Fye approach and have symmetrized the cluster 
Green function to obtain the corresponding quantity on 
the lattice: 

G(k,iu m )^ u = I^e*' (x - x ^G(K,ia; m ) M)M) . 

U a fi 8 

(2) 

Here x a denotes the unit cell positions within the super- 
cell, \i and v run over the c and / orbit als within each 
unit cell, and k and K differ by a reciprocal lattice vec- 
tor of the supercell Bravais lattice. The rotation to real 
frequencies was carrie d out with a stochastic analytical 
continuation technique! 12 * 13 ! 

III. RESULTS 

We consider the following model parameters: t c = tf = 
t, U c /t = Uf/t = 12, V/t = 1/2, e°Jt = 3, and e° f /t = 
0. We have chosen large values of U c and Uf to reflect 
the contact repulsion of the 3 He atoms and to guarantee 
that each single layer is well within the Mott insulating 
phase at half-band fillingpJ These values of the Hubbard 
interaction lead to low double occupancy, thus generating 
local moments. The difference e Q c — > is a crude 
accounting for the van der Waals forces (both 4 He- 3 He 
and 3 He- 3 He) that preferentially fill the first layer. 

A. Layer densities 

The generic Mott insulating state is characterized by 
a density (h) = 1 and a vanishing charge susceptibility: 
i-e., Xch = d(h)/d/jL = 0, where \i denotes the chemical 
potential. Figure [2] plots the layer-resolved densities (h c ) 
and (hf) as a function of the chemical potential, which 
controls the overall 3 He concentration. For both the odd 
and even cluster sizes, (hf) shows a plateau feature cen- 
tered around (hf) = 1, whereas (h c ) grows smoothly. In 
contrast to the Mott insulating state, xlh ~ d{hf)/dfi 
never vanishes. Hence, charge fluctuations between the 
layers are allowed and the simple picture of a com- 
plete decoupling of the layers never holds. Although 
the plateau feature is common to all cluster sizes, the 
data show distinct odd-even effects. For 7V U = 1 and 
N u = 3, (hf) is a continuous function of the chemical 
potential for all temperatures considered. In contrast, 
for the even clusters, N u = 2 and N u = 4, a disconti- 
nuity in (hf) emerges below a critical temperature T c 
and at a critical chemical potential. For N u = 4, a 
robust discontinuity is present at T c c± t/20, whereas 
for 7V U = 2 this feature already appears at T c ~ t/15. 
Since (h c + hf) = dF/dfi, where F is the free energy, 
the jump in the total fermionic density signals a density- 
driven first-order transition. In a canonical ensemble, 
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FIG. 2. The average occupation number is plotted for the upper- (top row) and lower-layer (bottom row) fermions as a function 
of chemical potential. Results are reported for both odd- (leftmost two columns) and even-numbered (rightmost two columns) 
clusters. Solid lines connect data sets with a common temperature. Vertical dashed lines mark the location of the emerging 
low-temperature discontinuity in (nf). 



states with total density lying within the jump are phase 
separated. 

B. First-layer effective mass and low temperature 
spectral functions 

We can estimate the / fermion's effective mass 
as a function of the chemical potential by con- 
sidering its cluster- averaged self-energy, E/(zo; m ) = 

Z)«=i s (/,a),(i>)(^m), and extracting the quantity 

t_ _ rn^ x _ _ Imjj^n = inT) 

T ~ m * ttT • U 

This estimate of the effective mass (or, equally, of the 
inverse of the coherence temperature To) is valid provided 
that the real space dependence of the self-energy is small 
and that the temperature T is extrapolated to zero. Data 
on the N u = 4 cluster presented in Ref . 1 1 shows that for 
fx < fJL c the assumption of a local self-energy is valid. 

At N u — 1 (see Fig.[3|, the effective mass increases as a 
function of chemical potential. This effect is also evident 
in the evolution of single particle spectral functions, 

A(k,oj) = -ImTrG(k,cj + zO + ), (4) 
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FIG. 3. The inverse of quantity defined in Eq. ([3]), which is 
proportional to the coherence temperature To, inverse effec- 
tive mass, and quasiparticle residue, is plotted as a function 
of chemical potential. 



plotted in Fig. |4l As exemplified by the data set at 
fi/t = -2.5 [FigTHTa)], the low -energy coherent features 
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FIG. 4. The amplitude of the single particle spectral func- 
tion, as defined by Eq. Q, is plotted for the N u = 1 cluster 
at various values of the chemical potential and temperature. 
The insets show the corresponding to slave-boson calculation. 
The inset to panel (a) shows the mean-field band structure 
consisting of two strongly-hybridized quasiparticle bands of 
mixed c and / character. In panel (b), the hybdridization 
is considerably weakened. In panel (c), the / fermions have 
dropped out of the Luttinger volume. 



of the spectral function compare favorably with a slave 
boson approximation leading to mass-renormalized hy- 
bridized bands. This state has a Luttinger volume that 
includes both / and c fermions, and the band with the 
largest Fermi volume has dominant / character. As a 
function of the chemical potential, the effective mass of 
the / band grows, and spectral weight is shifted to the 
upper Hubbard band. At \i — and fit = 30, the data 
of Fig. [4^b) exhibits typical heavy fermion character: 
a lower Hubbard band located at cjl/£ ~ —6, an up- 
per Hubbard band at ujjj jt — 6 = uj^jt + Uf /t, and a 
heavy band with dominant / character in close vicinity of 



the Fermi energy that hybridizes with a light conduction 
band. At N u = 3, the coherence temperature (see Fig.[3| 
is reduced with respect to the N u = 1 case but neverthe- 
less shows a similar overall behavior: a rapid decrease as 
a function of chemical potential followed by saturation 
at a lower value than for the N u = 1 case. Within the 
accessible temperature range of the N u = 3 cluster, the 
single particle spectral function shows the same features 
as for the N u = 1 case. 

For the even site clusters, the initial decrease of the 
coherence temperature is cut off by the first-order tran- 
sition. Far below /i c , as exemplified by fi/t = —2.5, the 
single particle spectral function is very similar to that 
observed on the N u = 1 cluster [cf. Figs. |4ja) and[5|a)]. 
With increasing chemical potential, the effective mass 
of the / band grows, and beyond /i c the / band drops 
out of the low-energy physics altogether. This can be 
understood at the static mean field level by a conven- 
tional slave boson theory in competion with local singlet 
formation in the first layer. The transition is signaled 
by the appearance of an anomalous expectation value 
Ay ~ (t}/U f ) EM&i,*)- The inset of Fig.gc) shows 
the band structure that results when this singlet order 
parameter breaks down the original lattice symmetry to 
that of the N u = 2 supercell. 

C. Spin susceptibilities and correlations 

We can extract from the cluster the local spin suscep- 
tibility as defined by 

X/(^ m ) = ^E /^e in ^<Sf(r).Sf(0)>. (5) 

A Fermi liquid below its Fermi temperature is Pauli para- 
magnetic, and hence Xfi^m — 0) is constant. On the 
other hand, a local moment is characterized by a Curie- 
Weiss law, Xf(i^m = 0) oc 1/(T + 0) at temperatures 
T ^> 0. As is apparent in Fig. [6j Xfi^m = 0) always 
exhibits a smooth crossover from the high-temperature 
Curie- Weiss to the low-temperature Pauli behavior, ir- 
respective of the lattice topology. For the odd lattice 
sizes, the crossover point tracks the coherence temper- 
ature. The same holds for the even lattice at \i < \i c . 
It is worth emphasizing that this qualitative change in 
magnetic response pinned to the coherence temperature 
(also denoted by To in Ref [7)) has been observed in the 
3 He bilayer experiment P Hence, at high temperatures, 
a local moment generated by the Hubbard interaction is 
present. The screening of this local moment, or in other 
words the quenching of its entropy, is at the origin of the 
different behavior between the odd- and even-numbered 
lattices. 

N u — 1. For this cluster smallest size, only the delocal- 
ized c fermions are available to screening the local mo- 
ment. This is precisely the Kondo effect, and one can 
view the heavy fermion paramagnetic state as originating 
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from the coherent, Bloch-like superposition of individ- 
ual Kondo screening clouds. Within a periodic Anderson 
model, this screening of the local moment is linked to a 
derealization of the / fermion. Hence, above To, when 
screening is absent, we expect the /-quasiparticle band 
to drop out of the low-energy physics. This is evident 
from Figs. [4^b) and |4^c) upon comparison of the high- 
and low-temperature spectral functions at \i = 0. 
N u = 3. This cluster size shows behavior very similar to 
that of the N u = 1 system, albeit with a lower coher- 
ence temperature. At values of the chemical potential 
where the / layer is approximately half filled, a mag- 
netic superexchange interaction J = 4t 2 /Uf is dynami- 
cally generated and the spin degrees of freedom on the 
first layer are described by a Heisenberg model on a three- 
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(a) A(k,u), fit = 20, fjb/t = -2.5 
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(b) A(k,o;), /3t = 20, fi/t = -1.875 
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FIG. 6. Temperature dependence of the inverse local spin 
susceptibility for the (a) N u = 3 and (b) iV u = 4 clusters. 
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FIG. 5. The amplitude of the single particle spectral function 
for the iV u = 4 cluster. The inset in panel (c) shows the single 
c-only band completely decoupled from the gapped, nearly 
flat band of the singlet-bound / fermions. The mean-field 
calculations were carried out on the N u = 2 cluster and the 
QMC on the N u = 4 systems. 



FIG. 7. (a) Normalized nearest-neighbor /-fermion spin-spin 
correlations for on the N u = 3 cluster. On a three-site ring, 
the value of this quantity for the Heisenberg model is given 
by —1/3. (b) Local dynamical spin structure factor. 
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site lattice. The ground state is fourfold degenerate cor- 
responding to a spin- 1/2 degree of freedom with either 
positive or negative chirality. As in the N u = 1 case, the 
only way to quench this residual entropy is via Kondo 
screening by the c fermions. To confirm this interpreta- 
tion, we have computed the nearest-neighbor spin-spin 
correlations on the first layer, S>(r) = (Sf • Sf +r ). Com- 
parison with the Heisenberg result is best achieved by 
normalizing the QMC data by the magnitude of the lo- 
cal moment, (r = 0). As shown in Fig. [7^ a), the 
nearest-neighbor antiferromagnetic spin-spin correlations 
are considerable. At high temperatures the energy scale 
at which they decay is set by the superexchange cou- 
pling J ft = 1/3; and at fi/t — —0.75, where we observe 
a Curie- Weiss law down to our lowest temperature, they 
compare favorably to the Heisenberg ground state result: 
S f (an)/S f (0) = -1/3. Fig. [^b) plots the dynamical 
local spin structure factor, 







S f (uj) =Im 



(6) 



at fit = 25 and as a function of the chemical potential. As 
mentioned above, at \ijt = —0.75 the residual entropy is 
not quenched. Consequently, a low frequency sharp fea- 
ture in Sf(u) marks the spin degenerate ground state of 
the three-site, half-filled Hubbard model. A feature at 
uj/t c± 0.5 corresponds to the first spin excitation, which 
for the three-site Hubbard model at Uf/t = 12 takes the 
value A sp /£ = 0.49. As we decrease the chemical po- 
tential from ji/t = —0.75 to ja/t = —1.5, the weight in 
the high-energy feature remains approximately constant, 
but the sharp low-energy feature decreases in intensity 
and is shifted to slightly higher energies. This screen- 
ing of the residual entropy by the c fermions competes 
with the nearest-neighbor antiferromagnetic fluctuations 
in the first layer and is at the origin of the upturn in 
S f (a.i)/S f (0) (see Fig. [7]) at low temperatures. 

Hence, at N u = 3 spin correlations between the / 
fermions quench part of the entropy associated with the 
formation of the local moments. The residual entropy is 
Kondo screened by the c fermions, and in comparison to 
the N u = 1 case leads to a suppressed coherence temper- 
ature. 

N u = 2, N u = 4 The even site clusters show a band- 
selective Mott transition and a low-energy decoupling of 
the first and second layers. As in the N u = 3 case, we 
can consider the effective Heisenberg model on the first 
layer. For even cluster sizes the ground state is unique 
and is spin singlet. The first-order transition we observed 
in Fig. [2] arises from competing screening mechanisms of 
the local moments generated by the nearly localized / 
fermions. On the one hand, the local moments can be 
Kondo screened by the light c fermions, thereby generat- 
ing heavy fermion behavior. On the other hand, they can 
form (among themselves) a spin singlet state entirely in 
the first layer. The gapping of the spin and charge degrees 
of freedom of the / quasiparticles at \i > \i c allows for 
the decoupling of / and c quasiparticles: a c quasiparticle 
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FIG. 8. Static and dynamical spin correlations on the iV u = 4 
cluster at f3t = 20. (a) Intracell spin-spin correlations between 
c and / fermions. (b) Nearest-neighbor spin-spin correlations 
between / fermions. (c) Local dynamical spin structure fac- 
tor. 



at the Fermi level cannot scatter off an / quasiparticle 
due to the absence of phase space. To support the pic- 
ture of a sudden change in the screening mechanism, we 
plot in Fig. [5] intra- and interlayer equal-time spin-spin 
correlations for the N u = 4 cluster size. At the critical 
chemical potential, we observe a sudden growth of the an- 
tiferromagnetic correlations between nearest-neighbor / 
fermions and a decrease in the intracell c-f spin-spin cor- 
relations. Figure [8] also shows the local dynamical spin 
structure factor. One observes a depletion of spectral 
weight at low energies on both sides of the transition and 
a considerable sharpening of the line shape in the band- 
selective Mott insulating state. At \i < fi c: we can inter- 
pret the data within an itinerant fermion picture where 
the mass enhancement prior to the band-selective Mott 
transition is taken into account by a renormalization of 
the hybridization V and hopping t as in a slave boson 
approachP Following this modeling, the peak position in 
Sf(uj) is expected to track the coherence temperature or, 
equivalently, the inverse effective mass. An explicit com- 
parison of those quantities is provided in Ref. [TTJ At 
H > /i c , Sf((j) should be interpreted within a localized 
/ fermion picture, in which case the peak position is a 
measure of the excitation energy required to break the 
singlet state of the / fermions. On a four-site Hubbard 
cluster, corresponding to the / layer in the N u = 4 case, 
this quantity is given by 0.214t and compares favorably 
to the data in Fig. [8] 



IV. CONCLUSIONS 

Our calculations are best summarized by the phase di- 
agrams plotted in Fig. [9J At high temperatures (T > T c ) 
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FIG. 9. Schematic phase diagrams for clusters that are (a) 
odd and (b) even in number. The solid line corresponds to the 
coherence temperature To, which sets the crossover scale be- 
tween states with large and small Fermi surfaces. The dashed 
line is a true phase boundary and corresponds to a line of 
first-order transitions terminating at a critical end point T c . 

the results are independent of the cluster size and topol- 
ogy. As the first layer approaches half band-filling the 
Hubbard U generates a large effective mass as well as 
local moments. The local moment is the key feature of 
the high temperature phase and results in a Curie- Weiss 
spin susceptibility. In the bilayer 3 He experiment this is 
indeed observed at temperatures above To prior to the 
solidification of the first layer. This high temperature 
local moment phase is characterized by a Luttinger vol- 
ume that counts only the c fermions. As the tempera- 
ture drops the entropy associated with the local moment 
has to be quenched, and competing quenching mecha- 
nisms are at the origin of the different phase diagrams. 
Heavy fermion character^ is associated with the screen- 
ing of the local moments by the conduction electrons. 
In the framework of the periodic Anderson model, the 
/ fermions delocalize so as to a generate the superex- 
change scale and reappear in the Luttinger sum rule. For 
our odd-sized cluster topologies only this scenario can oc- 
cur. It is important to note, however, that the step from 
N u = 1 to N u = 3 is linked to a substantial decrease 
of the coherence temperature since for the N u = 3 clus- 
ter the entropy is partially lifted due to intralayer spin 
correlations. 

On clusters of even size, the / fermions can form an 



insulating spin-singlet state and hence quench the en- 
tropy without involving the first layer fermions. This 
allows for a band-selective Mott transition — or Kondo 
breakdown — in which the / fermions drop out of the 
Luttinger volume down to the lowest temperature. De- 
spite the breaking of translation invariance inherent to 
the CDMFT, the Luttinger sum rule still holds when 
formulated in the Brillouin zone of the supercell Bravais 
lattice. 

Given this odd-even effect, the extrapolation to large 
cluster sizes is difficult and bound to be speculative. One 
can conjecture that for even site lattices, T c is set by the 
spin gap A sp of the corresponding half-filled Hubbard 
model of the first layer. At U f /t = 12, A sp /t = 0.325 for 
N u = 2 whereas A 8p /t = 0.214 for N u = 4. The decrease 
in T c between the N u — 2 and N u = 4 clusters is consis- 
tent with the decrease in the spin-gap. For odd lattices, 
one can follow the idea that the coherence temperature 
tracks the residual entropy per site of the half-filled Hub- 
bard model on the first layer. Given the above conjecture 
and the fact that the Hubbard model on a triangular lat- 
tice has a unique ground state, we arrive at the conclusion 
that the coherence temperature indeed vanishes beyond 
a critical chemical potential. This stands in agreement 
with the slave boson calculations of Ref. [8] If the mag- 
netic system on the first layer orders and breaks a lattice 
symmetry, then the /-fermions can drop out of the Lut- 
tinger volume without violating the Luttinger theorem. 
The stability of such a phase with respect to a finite hy- 
bridization matrix element generating a Kondo coupling 
between the layers has been discussed in Ref. [T6l On the 
other hand, if no symmetries are broken such that a spin 
liquid state is realized on the first layer, 17 fractionalized 
Fermi liquids as proposed in Refs. [18] and [19] could be 
realized. 

It is interesting to recast our results in terms of the Q- 
K phase diagram for heavy fermions proposed by Cole- 
man and Nevidomskyy! 2 ^ Here, K corresponds to the 
magnitude of the Kondo screening and Q is a measure 
of frustration between the / fermions. For our model in 
the local moment regime, the Kondo coupling between 
the two layers is dynamically generated starting from 
second-order perturbation theory in the hybridization. 
The frustration between the / fermions is generated by 
the hopping matrix element tf, which again in the lo- 
cal moment regime leads to a superexchange interaction 
between the / fermions. In the framework of Ref. [20l 
and in agreement with our numerical simulations, it is 
the frustration between the / fermions that drives the 
Kondo breakdown or band-selective Mott transition. For 
a recent review in this domain, we refere the reader to 
Ref. [21] This is in contrast to the Kondo lattice model, 
where the / fermions interact solely through the RKKY 
interaction, an d no K ondo breakdown is observed in clus- 
ter simulations! 22 * 23 * 

Let us finally return to the bilayer 3 He experiment. 
Heavy fermion character is clearly seen by the increase 
of the effective mass (or decrease of the coherence tern- 
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perature). Furthermore, and as seen in our calculations, 
To marks the crossover between a Curie- Weiss and Pauli 
behavior of the spin susceptibility. The Q or compet- 
ing interaction which localizes the / fermions to the first 
layer are the three-body exchange processes. These pro- 
cesses, which in solid 3 He can dominate the two-body an- 
tiferromagnetic exchange, lead to the observed ferromag- 
netic behavior. Given this interpretation of the experi- 
ment, an extremely important issue would be to pin down 
the experimental value of the Weiss constant. Above 
T it should be positive and essentially track the Kondo 
scale. The transition to the ferromagnetic state should 
be accompanied by a vanishing and subsequently nega- 



tive value of the Weiss constant. 
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